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"It  is  out  of  the  question  to  formulate  and  carry  out 
experiments  nowadays  without  leaking  use  of  similarity 
and  dimensionality  concepts.11 


I . INTRODUCTION 

One  measure  of  the  importance  of  a scientific  area  is  the  number 
of  papers  written9in  that  area.  Dimensional  analysis  has  generated  over 
600  publications,  however,  only  a handful  of  these  papers  deal  with 

computerized  dimensionality . ^ 6, 7 

There  are  several  reasons  why  so  little  computer  work  has  been  done 
in  dimensional  analysis.  First,  much  work  has  been  devoted  to  method- 
ologies and  specific  solutions  in  highly  specialized  areas.  Once  this 
work  has  been  done,  a computer  is  not  necessary.  Second  and  more  impor- 
tant is  that  computer  codes  .have  been  inflexible.  Only  one  solution  is 
usually  generated  for  a given  ordering  of  dimensional  equations.  There- 
fore, fruitful  and/or  well  recognized  solutions  might  not  be  generated. 
This,  just  because  of  ordering.  Next,  identifiable  and  well  recognized 
solutions  are  not  identified.  The  user  does  needless  computational  work 
when  this  occurs.  Finally,  with  the  exception  of  the  work  of  Cohen  and 
Ferrari  (Reference  6),  solutions  are  not  usually  numerically  evaluated. 
Therefore,  a great  deal  of  computational  work  must  be  performed  before 
suitable  data  is  available  for  plotting,  regression,  clustering,  or 
factor  analysis. 


^Sedov,  L.  I.  Similarity  and  Dimensional  Methods  in  Mechanics.  New  York, 
Academic  Press,  1959. 

2 

“Happ,  W.  W.,  Private  Communication,  January  1976. 

3Happ,  W.  W.,  "Computer-Oriented  Procedures  for  Dimensional  Analysis." 
Journal  of  Applied  Physics,  38,  3918-3926,  September  1967. 

4Sloan,  A.  D.  and  Happ,  W.  W.,  "Computer  Program  for  Dimensional  Analysis.' 
Electronics  Research  Center,  Cambridge,  Mass.,  NASA  TN  D-5165,  April  1969 

°Chen,  W.  K.  "Algebric  Theory  of  Dimensional  Analysis."  Journal  of  the 
Franklin  Institute,  292,  403-422,  December  1971. 

6Cohen,  J.  and  Ferrari,  J.  0.,  "A  Conversational  Language  for  Solving 
Problems  in  Dimensional  Analysis."  Computer  Methods  in  Applied  Mechanics 
and  Engineering,  5,  53-67,  January  1975. 

7Hirschberg,  M.A.,  "A  Computer  Solution  of  the  Buckingham  Pi  Theorem 
Using  SYMBOLANG,  A Symbol  Manipulation  Language,"  USA  Ballistic  Research 
Laboratories  Report  No.  1824,  August  1975.  (AD  #A016901) 
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The  program  described  in  this  paper  overcomes  most  of  the  named 
deficiencies > that  is,  those  of  inflexibility.  Computer  solutions  are 
generated  using  the  Buckingham  Pi  Theorem  (Section  II).  Computations 
are  performing  using  SYMBOLANG  (Section  III),  an  algebraic  (symbol)  mani- 
pulation package.  Solutions  to  dimensional  equations  are  kept  in  symbo- 
lic form.  In  addition,  solutions  are  manipulated  in  two  ways.  First, 
many  allowable  combinations  of  pi  terms  (up  to  720)  can  be  generated  by 
the  program  automatically  changing  the  order  of  the  dimensional  equations 
Second,  a small  group  of  operators  (squaring,  cubing,  square-rooting, 
cube-rooting,  reciprocating,  and  division  by  both  the  minimum  and  maxi- 
mum exponent)  are  applied  to  each  generated  solution.  Various  transfor- 
mations applied  to  solutions  produce  new  solutions.  Products  of  solutions 
also  produce  new  solutions.  Provision  has  been  made  for  numerical  sub- 
stitution into  the  symbolic  forms  producing  scalars  (evaluated  solutions) 
for  subsequent  use.  In  addition,  every  symbolic  solution  is  compared 
with  well-known  nondimensional  numbers  (e.g.,  Reynold's,  Weber's,  etc.). 
Well-known  numbers  are  selected  from  the  Land  Table  . When  a match  be- 
tween a symbolic  and  a Land  number  is  found,  the  name  of  the  Land  number 
is  printed  for  all  such  numbers  which  match  (i.e.,  there  may  be  more 
than  one  name  printed  for  a particular  solution) . 

II.  THE  BUCKINGHAM  PI  THEOREM 
9 

The  Buckingljigm  Pi  Theorem  summarizes  the  entire  theory  of  dimen- 
sional analysis.  The  result  of  a dimensional  analysis  is  the  reduction 
of  the  number  of  variables  in  a problem.  Application  of  the  Pi  theorem 
itself  provides  the  method  of  solution  of  a set  of  dimensional  equations. 
Simply  stated,  the  pi  theorem  asserts: 

If  there  are  n variables  involving  N fundamental  units, 
these  may  be  combined  to  form  n-N  dimensionless  parameters 
each  involving  N + 1 variables. 

The  usual  method  of  applying  the  Pi  theorem  is  for  one  to  write  equa 
tions  describing  the  physical  system  one  is  interested  in  in  terms  of  a 
set  of  fundamental  units  (force,  length,  time,  angle,  or  mass,  length 
time,  etc.).  The  equations  are  then  systematically  exponentiated  and 
multiplied  together  (hence  "Pi"  theorem  for  the  mathematical  symbol  for 
multiplication  (tt)  . The  resulting  exponentials  form  an  N by  N set  of 


g 

Land,  N.  S.,  "A  Compilation  of  Nondimensional  Numbers,"  Washington,  DC, 
US  Government  Printing  Office,  NASA  SP-274,  1972. 

9 

Buckingham,  E.,  "On  Physically  Similar  Systems:  Illustrations  of  the  Use 
of  Dimensional  Equations,"  Physical  Review,  2_,  345-376,  1914. 

^Langhaar,  H.  L.,  Dimensional  Analysis  and  Theory  of  Models,  New  York, 
Wiley,  1951. 
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linear  equations  whose  solution  is  applied  back  to  the  variables  of  the 
problem  (see  the  example  below) . One  can  see  that  it  is  tedious  to  work 
a problem  involving  many  variables  by  hand.  In  fact,  there  is  no  good 
indicator  of  just  how  many  variables  one  should  include  in  a problem 
(this  is  true  even  for  a computerized  solution) . 

The  following  example  (taken  from  Housner  and  Hudson11)  used 
throughout  the  remainder  of  the  paper  illustrates  the  mechanisms  of  the 
Buckingham  Pi  Theorem. 

Consider  a drag  force  (F)  acting  upon  a body  moving  through 
a fluid.  Assume  a constant  velocity  (v)  through  the  fluid 
of  density  (p)  and  viscosity  (y) . If  the  analysis  is  applied 
to  bodies  of  a specific  shape,  the  cross-sectional  area  (A) 
may  be  used  as  a measure  of  the  body's  size. 

The  following  variables  and  fundamental  units  enter  into  the  problem: 


Variable 

Fundamental  Units 

F 

= F 

-2m 

U 

= FL  T 

2 

A 

= L 

-4  2 

P 

= FL  T 

-1 

V 

= LT 

where  F = Force,  L = Length  and  T = Time, 

According  to  the  Pi  theorem,  two  terms  can  be  formed  from  the  five  equa- 
tions (each  equation  is  expressed  in  terms  of  the  three  fundamental  units) . 
The  two  solutions  will  each  contain  four  of  the  variables.  The  pi  terms 
formed  with  this  ordering  are: 


ir  = FA°0Vy  = 


pl  + eL2  0f-4B+YT2g-y 


and 


WA  “V  = 


pl+3L~2+2  a-43+Y^l+23-Y 


Solving  these  equations  we  find:  has  the  solution  a=  -1,  0 = -1, 
Y = -2;  u-  has  the  solution  a=  -1/2,  0 = -1,  y = -1;  so  the  result- 
ing dimensionless  pi  terms  are: 

^Housner,  C.  W.  and  Hudson,  D.  E.,  Applied  Mechanics  Dynamics.  New 
York,  van  Nostrand,  1950. 
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F 


and 


y 


TT,  = 


ApV 


A%V 


7T  is  a pressure  coefficient  and  tt9  is  the  reciprocal  of  the  Reynold's 

i 10  ^ 

number  . One  can  see  how  the  equations  are  systematically  selected  and 
visualize  the  ease  with  which  such  an  algorithm  can  be  programmed. 

The  above  represent  one  set  of  solutions.  In  this  particular  case, 
eight  other  solutions  are  possible  depending  upon  the  reordering  of  the 
five  basic  equations.  The  full  set  of  solutions  is  listed  below: 


(1) 


J? 

ApV2 


(6) 


ApV2 

F 


(2) 


(3) 


2 2 
Ay  V 


(8) 


y (A°) 


pF  CV°) 

2 


(9) 


VyA2 

F 


(5) 


y (V°) 


(10) 


pF(A°) 

2 

y 


In  this  example,  only  solutions  (1)  and  (2)  are  fruitful  (contain 
information) . The  other  eight  solutions  can  all  be  derived  from  solu- 
tions (1)  and  (2) . We  must  note,  however,  that  solutions  (4)  and  (10) 
are  minimal  in  the  sense  that  they  (a)  contain  the  fewest  number  vari- 
ables and  (b)  the  sum  of  their  (integer)  exponents  is  a minimum  (see 
References  4 and  5 for  a further  discussion  of  minimal  solutions) . 
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III.  SYMBOLANG 


SYMBOLANG1  5 , a high-level  FORTRAN  language  for  algebraic  (symbol) 

manipulation,  was  used  to  form  the  symbolic  products  (pi  terms)  of  the 
fundamental  units.  This  application  was  well  suited  for  SYMBOLANG.  Not 
only  were  solutions  generated,  but  one  was  able  to  see  the  development 
of  the  pi  terms  as  the  products  were  being  formed.  In  addition,  the 
final  equation  was  also  displayed.  SYMBOLANG  is  of  pedagogical  value 
in  demonstrating  the  mechanisms  of  the  Buckingham  Pi  Theorem. 


IV.  THE  PROGRAM 


Listings  of  the  routines*  described  in  this  section  appear  in  Appen- 
dix I.  The  many  SLIP*  and  SYMBOLANG*  routines  are  not  included;  however, 
they  are  available  from  the  authors  of  Reference  12. 

The  name  of  the  main  program  is  BUCKY.  BUCKY  initializes  the  SLIP- 
SYMBOLANG  working  storage  area,  then  reads  and  displays  the  inputs 
(sample  inputs  appear  in  Appendix  II  while  an  explanation  of  the  inputs 
and  formats  required  for  them  appears  in  Appendix  III).  BUCKY  next  cal- 
culates the  number  of  combinations  (orderings)  possible  based  upon  the 
number  of  equations  input;  however,  only  720  permutations  are  allowed. 
BUCKY  next  breaks  the  equations  into  two  pieces;  the  ones  to  the  left  of 
the  equal  sign  are  variables  and  those  to  the  right  of  the  equal  sign  are 
fundamental  units.  Now,  the  Pi  theorem  algorithm  is  invoked.  This  is 
where  the  exponentiation  and  multiplication  of  fundamental  units  occur. 
The  result  of  this  step  is  the  system  of  equations  which  is  solved 
using  a matrix  inverse  routine.  The  solution  of  the  linear  equations 
is  performed  in  subroutine  BUCKSV.  Next,  the  appropriate  numerical 
operator  (squaring,  etc.)  is  established  and  a numerical  substitution 
made.  A scalar  is  produced. 

* Variable  length  calls  and  calls  to  functions  have  been  modified  to  run 
on  the  BRLESC  II  Computer. 


12 

Findler,  N.  V.,  Pfaltz,  J.  L.,  and  Bernstein,  H.  J.,  Four  High-Level 
Extensions  of  FORTRAN  IV;  SLIP,  AMPPL-II,  TREETRAN,  SYMBOLANG.  New 
York,  Spartan,  305-387,1972. 

1 3 

Hirschberg,  M.  A.,  nSYMBOLANG  - A SLIP  Extension  for  Algebraic  Manipu- 
lation," USA  Ballistic  Research  Laboratory  Report  No.  1749,  November 
1974.  (AD  #A003190) 


9 


Evaluation  is  performed  in  subroutine  EVALP.  EVALP  is  an  inelegant 
routine  whose  main  virtue  lies  in  producing  the  correct  scalar  value.  Up 
to  five  values  may  be  input  for  each  of  24  variables  (a  considerable  num- 
ber for  a dimensional  analysis) . Reciprocating  is  also  performed  in 
EVALP . 

After  evaluation,  the  exponents  of  the  symbolic  solution  are  passed 
to  TABLUK,  which  compares  them  with  well-known  dimensionless  numbers 
taken  from  the  Land  Table.  When  a match  is  found,  the  name(s)  of  all 
numbers  in  the  table  fitting  the  match  are  printed. 

Following  the  table  lookup,  a new  numerical  operator  is  selected 
(subroutine  PRMTE)  and  the  evaluation  and  lookup  process  repeated.  When 
all  numerical  operators  have  been  processed,  a new  ordering  of  the 
equations  is  generated  (this  is  a random  process  performed  in  ONEMR 
which  also  keeps  track  of  which  permutations  have  already  been  made) 
and  the  entire  solution  process  is  repeated.  Up  to  720  solutions  are 
permitted,  so  all  solutions  will  be  generated  even  for  relatively  large 
problems.  Sample  program  outputs  appear  in  Appendix  IV. 


V.  DISCUSSION 


In  an  earlier  section,  we  have  seen  what  program  BUCKY  does.  Namely, 
it  forms  solutions  to  the  Buckingham  Pi  Theorem,  evaluates  these  solu- 
tions, and  finally,  identifies  them  when  possible.  In  addition,  reorder- 
ing the  dimensional  equations  allows  for  different  solution  sets  to  be 
formed.  When  there  are  few  equations  to  be  solved,  every  possible  solu- 

3 4,5 

tion  is  formed;  therefore,  an  optimal  solution  is  always  generated  , ' . 
Furthermore,  numeric  substitution  into  the  symbolic  form  generates  scalar 
values  (evaluated  solutions)  which  may  be  kept  in  a data  base  for  further 
use.  When  an  investigator  finds  a solution  particularly  suited  to  his 
needs,  selective  retrieval  of  evaluated  solutions  allows  data  to  be  plot- 
ted as  well  as  used  in  regression,  clustering,  or  factor  analysis.  As 
a last  step,  an  attempt  is  made  to  identify  each  solution  (primary  and 
algebraically  manipulated)  by  comparing  solutions  with  well-known  non- 
dimensional  numbers.  It  should  also  be  noted  that  by  keeping  solutions 
in  symbolic  form  one  does  not  need  to  refer  to  other  documents  to  deter- 
mine which  variables  are  involved  in  the  solution. 

There  are  some  shortcomings  to  program  BUCKY.  First,  no  attempt 
is  made  to  verify  that  there  is  a consistent  set  of  equations  (none  of 
the  other  computer  programs  does  this  either) . Failure  to  provide  a 
consistent  set  of  equations  results  in  a singular  matrix  which  cannot  be 
inverted.  Second,  the  program  does  not  provide  for  an  automatic  change 

of  units6.  This  is  correctable  by  providing  conversion  factors  and 
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having  a separate  conversion  calculation  phase  before  evaluation  occurs. 
Finally,  while  no  other  program  attempts  to  identify  solutions,  not  all 
well-known  solutions  are  tabled.  This  can  be  easily  corrected  by  add- 
ing solutions  to  the  table.  In  addition,  some  well-known  solutions  are 
not  recognized  because  of  dimensional  substitutions  which  can  be  but  are 
not  made.  For  instance, 

A 2pV 

y 


is  a Reynold's  number;  however,  the  program  does  not  realize  that  A = L 
(A  is  area,  L is  length) . This  fault  can  also  be  corrected  by  substitut- 
ing fundamental  units  for  variables;  however,  the  expense  of  making  all 
such  substitutions  does  not  seem  to  warrant  the  gain  (there  are  always 
numerous  little  things  one  can  do I) 

No  further  extension  of  this  work  is  contemplated  at  the  present 

time . 


ACKNOWLEDGEMENT 


The  author  would  like  to  thank  Dr.  Benajmin  E.  Cummings  for  encour- 
aging this  work  and  for  providing  a meaningful  data  set  for  the  example 
used  in  this  report. 


11 


REFERENCES 


1.  Sedov,  L.  I.,  Similarity  and  Dimensional  Methods  in  Mechanics,  New 
York,  Academic  Press,  1959. 

2.  Happ,  W.  W.,  Private  Communication,  January  1976. 

3.  Happ,  W.  W.,  "Computer-Oriented  Procedures  for  Dimensional  Analysis", 
Journal  of  Applied  Physics,  38,  3918-3926,  September  1967. 

4.  Sloan,  A.  D.  and  Happ,  W.  W.,  "Computer  Program  for  Dimensional 
Analysis",  Electronics  Research  Center,  Cambridge,  MA,  NASA  TN 
D-5165,  April  1969. 

5.  Chen,  W.  K.,  "Algebric  Theory  of  Dimensional  Analysis",  Journal  of 
the  Franklin  Institute,  292 , 403-422,  December  1971. 


6.  Cohen,  J.  and  Ferrari,  J.  0.,  "A  Conversational  Language  for  Solv- 
ing Problems  in  Dimensional  Analysis",  Computer  Methods  in  Applied 
Mechanics  and  Engineering,  5_,  53-67,  January  1975. 

7.  Hirschberg,  M.  A.,  "A  Computer  Solution  of  the  Buckingham  Pi  Theorem 
Using  SYMBOLANG,  A Symbol  Manipulation  Language",  USA  Ballistic 
Research  Laboratory  Report  No.  1824,  August  1975.  (AD  #A016901) 

8.  Land,  N.  S.,  "A  Compilation  of  Nondimensional  Numbers",  Washington, 
DC,  US  Government  Printing  Office,  NASA  SP-274,  1972. 

9.  Buckingham,  E.,  "On  Physically  Similar  Systems:  Illustrations  of 
the  Use  of  Dimensional  Equations",  Physical  Review,  2_,  345-376,  1914. 

10.  Langhaar,  H.  L.,  Dimensional  Analysis  and  Theory  of  Models,  New  York, 
Wiley,  1951. 

11.  Housner,  C.  W.  and  Hudson,  D.  E.,  Applied  Mechanics  Dynamics,  New 
York,  van  Nostrand,  1950. 

12.  Findler,  N.  V.,  Pfaltz,  J.  L.  and  Bernstein,  H.  J.,  Four  High-Level 
Extensions  of  FORTRAN  IV,  SLIP,  AMPPL-II,  TREETRAN,  SYMBOLANG.  New 

York,  Spartan,  305-387,  1972. 

13.  Hirschberg,  M.  A.,  "SYMBOLANG,  A SLIP  Extension  for  Algebraic 
Manipulation" , USA  Ballistic  Research  Laboratory  Report  No.  1749, 
November  1974.  (AD  #A003190) 


12 


APPENDIX  I 

PROGRAM  COMPUTER  LISTING 


13 


PROGRAM  BUCKY 


C 

C THIS  IS  THE  PI  THEOREM  SOLVER 
C WRITTEN  BY  MA  HIRSCHBERG 
C oCTCBER  1974 
C 

C SET  UP  SLIP  STORAGE 

common  avsl,x<ioo) 

COMMON/MAXMIN/XAMAX,XAMIN 
DIMENSION  SP ( 20000 ) 

C SET  UP  PROGRAM  STORAGE 

DIMENSION  LNEW(  100  ) , LTEMP ( ICC ) ,LPGWER(100) 

DIMENSION  L FORMS ( 100) 

DIMENSION  APRIMSI 6, 100) 

DIMENSION  XFRMLA( 100 ) ,XSOLS( ICO) 

DIMENSION  I V AL ( 24 ) 

DIMENSION  SOLTN ( 24 ) 

DIMENSION  I Z EE ( 78 ) 

DIMENSION  FFCRM(IOU) 

DIMENSION  L EXP ( 78 ) 

DIMENSION  EEXPO ( 5 ) 

DIMENSION  NSGLS(720) ,IRAYS(20) 

DATA  EEXPO/ lo , 2»  , 3. «• 5, « 3333333333/ 

DATA  IMAXY/7 / 

C DEFINE  EXPONENTS  A-Z»A1-Z1,A2-Z2 

DATA  IZEE/1HA, 1HB.1HC, 1HD , 1HE , 1HF , 1HG , 1HH , 1H I , 1H J , 1HK , 1 HL , 1HM, 1HN, 

1 1H0* 1HP , 1HQ , 1HR , lHStlHT  » 1HU, 1HV, 1HW, 1HX, 1HY,1HZ, 

22HA1,2HB1,2HC1,2HD1,2HE1,2HF1,2HG1,2HH1,2HI1,2HJ1,2HK1, 2HL1,2HM1, 
32HN1,2HJ1,2HP1,2HQ1,2HR1,2HS1,2HT1,2HU1,2HV1 ,2HW1 ,2HX1 , 2HY1 , 2HZ1 , 
42HA2,2HB2,2HC2, 2HD2, 2HE2 , 2HF2 » 2HG2 » 2HH2 , 2H I 2 , 2H J2 , 2HK2 , 2HL2 , 2HM2 , 
52HN2,2H02,2HP2, 2HQ2, 2HR2, 2HS2 , 2HT2 , 2HU2 , 2HV2 , 2HW2 ,2HX2 , 2HY2 , 2HZ2/ 

C 

C SETUP  WORKING  STORAGE 

CALL  INITASI SP, 20000) 

2 CONTINUE 

C READ  NUMBER  OF  PRIMITIVES 
READ  (5,10)  NPRMS 
10  FORMAT  (15) 

WRITE  (6,12)  NPRMS 

12  FORMAT( 1H1, 55X, 18H  PI  THEOREM  SOLVER/ 

1 2 1H  NUMBER  OF  PR IM I TI VE S , 4X , I 5 ) 

WRITE  (6,13) 

13  FORMAT  (26H0LAND  CANDIDATES  REFERS  TO/ 

1 57H  A COMPILATION  OF  DIMENSIONLESS  NUMBERS  BY  NORMAN  S.  LAND/ 

2 48H  NASA  SP-274,  USGPO  STOCK  NUMBER  3300-0408,  1972) 

C THE  PRIMITIVES  ALSO  CONTAIN  VALUES  FOR  DETERMINING  THE  VALUE  OF  A PI 
C TERM;,  IF  A VALUE  IS  -9999.9,  THEN  NO  MORE  EVALUATIONS  WILL  BE  MADE  WITH 
C THAT  VARIABLE. 

WRITE  (6,14) 

14  FORMAT  ( 30HUPR I M I T I V ES  USED  IN  EVALUATION) 

WRITE  (6,18) 

18  FORMAT  ( 1H  , 10HPR  IM  I T I VE S , 22X , 6HVALUES ) 

C READ  PRIMITIVES 
IE VAL=u 

DO  25  1=1, NPRMS 

RE  AC  (5,21)  (APRIMS(J, I ) , J= 1 , 6 ) 
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21  FORMAT  ( A10, 2X,  5E12.8) 


COUNT  THE  NUMBER  OF  VARIABLES  FOR  EVALUATION 
11=0 

DO  22  J = 2 * 6 

IF  ( APRIMS( J, I ).EQ.-9999.9)  GO  TO  23 

11=11+1 

IEVAL=1 

22  CONTINUE 

23  CONTINUE 
IV  AL ( I ) = 1 1 

OUTPUT  PRIMITIVES 

WRITE  (6,24)  ( APRIMS( J, I ) ,J=1,6) 

24  FORMAT  ( 1H  , A 10 , 2X, 5E20. 8 ) 

25  CONTINUE 

READ  NUMBER  OF  FORMULAS 
READ  (5,10)  NFRMS 

READ  NUMBER  OF  INDEPENDENT  VARIABLES 
READ  (5,10)  IVARS 
WRITE  (6,32)  NFRMS, IVARS 
32  FORMAT  ( 10HOTHERE  ARE, 1X,I5,1X, 8HF0RMUL  AS , 

1 1X,9H INVOLVING, IX, 15, IX , 9HVAR I ABLES  ) 

I X TRA  = I V ARS  + 1 

CALL  COMBY(NFRMS, IXTRA,ICOMS) 

WRITE  (6,34)  ICOMS 

34  FORMAT  ( 10H0THERE  ARE,  IX,  I 5,  IX,  211+POSSIfiLE  COMBINATIONS) 
ICCC=0 

DO  36  1=1, NFRMS 
IR  AYS ( I ) = I 
ICCC=10*ICCC+I 
36  CONTINUE 

NSOLS ( 1 ) = ICCC 
NCT=  1 

READ  FORMULAS 

DO  40  1=1, NFRMS 

CALL  INL  1ST ( L FORM S( I), 5H INPUT, 3HVAL, TEMP) 

40  CONTINUE 

CALCULATE  THE  NUMBER  OF  PI  TERMS 
IP  I = NFRMS-IVARS 
WRITE  (6,42)  IPI 

42  FORMAT  ( 10H0THERE  ARE, IX, I 5, IX, 8HPI  TERMS) 

OUTPUT  FORMULAS 

DO  50  1=1, NFRMS 

CALL  LSQPNT ( LFORMS ( I ) , 8HF0RMULAS, 999. , TEMP ) 

5u  CONTINUE 

STRIP  EQUALS  OFF  AND  SETUP  NEW  SYMBOLANG  LISTS 
IG0N=0 
55  CONTINUE 
IX  DT  = 0 
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DO  200  1=1* NFRMS 

SET  UP  READER  FOR  FORMULAS 
II I=IRAYS( I ) 

LRD=LRDR0V( L FORMS (III)) 

STRIP  FIRST  PART  OF  FORMULA 
DO  60  1 1 = 1*  A 

DATUM=ADVSER( LRD, FLAG) 

IF  (I  Io  N E*  3 ) GO  TO  60 
IXDT=IXDT+1 

SAVE  VARIABLE  NAME  IN  FORMULA  FOR  LATER  USE 
FFORM ( IXDT ) =CATUM 
60  CONTINUE 
LW=0 

LW=LIST( L W ) 

IC  = 0 

63  CONTINUE 

ACVANCE  THROUGH  LIST 

DATUM=ADVSER( LRD, FLAG) 

I F ( FLAG.NE.O. ) GO  TO  70 

SET  UP  TEMPORARY  LIST  AND  COUNT  ELEMENTS 
IC=IC+1 

CALL  NEWBOT( DATUM, LW, TEMP ) 

GO  TO  65 
7v  CONTINUE 

SET  COUNT  AND  FORM  NEW  NEW  SYMBCLANG  LIST 
LN  EW ( I ) =0 

LN  EW ( I )=LIST(LNEW( I ) ) 

LRD=LRDRuV( LW ) 

LT  EMP ( I ) =L I ST  C 9 ) 

CALL  NEWBOT(LTEMP( I ) ,LNEW(  I ),TEMP) 

LC= I C-2 

IF  (LCoLEoO)  GO  TO  2C00 
DO  80  J=  1 » LC 
DATUM=ADVSER(LRD, FLAG) 

CALL  NEWBOT( DATUM, LTEMP( I ) ,TEMP) 

80  CONTINUE 

ERASE  TEMPORARY  LIST  AND  PRINT  NEW  LIST 
CALL  I RALST ( LW , TEMP ) 

C CALL  LSQPNT(LNEW( I ) , 4HLNEW, 999. , TEMP ) 

20 U CONTINUE 

IF  ( IGON.EQ.l)  GO  TO  220 
I GON= 1 
K=NFRMS- 1 


(SOLUTION) 
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DO  210  1=1, K 


SETUP  EXPONENTS 

LEXP ( I ) = LSQMN1( LSQMN3( 1* , IZEE ( I ) , 1.  ) ) 

2 1C  CONTINUE 
220  CONTINUE 

SETUP  PI  TERMS 

LFRMS=LSCMN 1 ( LSQMNK 1.  ) ) 

DO  260  1 = 1, IPI 
LTV=LSQMNll  LSQMNK  1*  ) ) 

KK=0 

CYCLE  THROUGH  FORMULAS  MOST  IMPORTANT  AND  SKIP  FORMULAS  SO  AS 
TO  INCLUDE  ONLY  THE  IMPORTANT  FORMULAS  ONE  AT  A TIME 
DO  245  J=1,NFRMS 
IF  (J.EQ.I)  GO  TO  215 
IF  ( J+IVARS.LE.NFRMS)  GO  TO  245 
KK>=KK  + 1 


RAISE  POWER 

LPOWER(KK)=LSQRAZ(LNEW(J) ,LEXP(KK) ) 

LTU  = LSQMEX(LTV,LPOWER(KK)  ) 

CALL  LSQDES ( LTV, TEMP ) 

TLGD=SEQRDR(LTU) 

LT V=LSQCPY (TLGD ) 

CALL  LSQDES( LTU,TEMP  ) 

GO  TO  245 
215  CONTINUE 

CALL  LSQDES( LTV, TEMP  ) 

LTV=LSQMEX(LFRMS,LNEW(J) ) 

245  CONTINUE 

CALL  LSQPNT( LTV , 3HLTV, 99^. ,TEMP) 

CALL  PRLSTS ( LTV , 4 ) 

ScLVE  EQUATIONS 

CALL  BUCKSV (LTV,IZEE,SOLTN,ICCUNT) 

CALL  LSQDES( LTV, TEMP ) 

DO  250  J = 1 , KK 

CALL  LSQDE$(LPOWER(J ) ,TEMP) 

250  CONTINUE 

PRINT  SOLUTION  FOR  PI  TERM 
IT IM  E = 1 

251  CONTINUE 
N0M=O 

WRITE  (6,252) 

252  FORMAT  ( 45H0 SOLUT I ON  OR  MANIPULATED  SOLUTION  FOR  PI  TERM) 
KK  =0 
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DO  258  J = l,  NFRMS 

IP  (J.NE.I)  GO  TO  255 

IF  (ITIMc.LE.5)  PTOUT=EEXPO ( I TIME ) 

IP  (ITIME.EQ.6)  PTOUT=l./XAMIN 
IP  (ITIME.EQo7)  PTOUT=l. /XAMAX 
WRITE  (6,253)  FFORM ( I ) , PTOUT 
NOM=NOM+l 

XPRMLA(NOM)=FFORM( I ) 

XSOLS(NOM)=PTOUT 

253  FORMAT  ( IH  , A 10, IX, 2H**, IX, El 4.8) 

255  CONTINUE 

IF  (J+IVARS.LE. NFRMS)  GO  TO  258 
KK=KK+ I 

WRITE  (6,253)  FFORM ( J ), SOLTN ( KK ) 

NO  M = NOM+  1 

XFRMLA(NOM)=FFORM( J) 

XSOLS(NOM)=SCLTN(KK) 

258  CONTINUE 

IP  (ITIME.GE.6  .AND.  I FI  X ( ABS ( PTOUT) ) • EQ. 1 ) GO  TO  259 
IF  (IEVAL.EQol) 

1CALL  EVALP( APR  I MS , XFRMLA , XSOL S, NOM, NPRMS, I VAL ) 

CALL  TABLUK ( NOM , XSOL  S ) 

259  CONTINUE 

IF  ( ITIME.GEo IMAXY)  GO  TO  260 
CALL  PRMTE( ITIME, IVARS, SOLTN) 

IT  IME  = IT IME+ I 
GO  To  251 

260  CONTINUE 

DO  270  1=1, NFRMS 

CALL  IRALST ( LTEMP ( I ) ,TEMP) 

CALL  IRALST(LNEW(  I), TEMP) 

270  CONTINUE 

CALL  ONEMR ( NSOL S , IRAYS, NCT, NFRMS, ICOMS) 

IF  (NCToGE.l)  GO  TO  55 
C Gr  TO  2 

CALL  EXIT 
2000  CONTINUE 

WRITE  (6,2010) 

2010  FORMAT  ( 7H  NO  NO  ) 

CALL  EXIT 
END 


19 


oo  oo  ooooo 


SUBROUTINE  BUCKSVtLIST,  ZEE , SOLTN, ICOUNT) 

THIS  SUBROUTINE  SETS  UP  THE  SOLUTION  FOR  THE  PI  THEOREM 
WRITTEN  BY  MA  HIRSCHBERG 
JANUARY  1975 

DIMENSION  Z EE ( 78  > 

DIMENSION  SOLTN ( 24 ) 

DIMENSION  AM AT (24,24) 

DIMENSION  AWORD ( 3 ) 

CLEAR  STORAGE 
DO  5 1=1,24 
S0LTN(I>=0.0 
DO  4 J=1 , 24 
AMAT( I, J)»0.0 

4 CONTINUE 

5 CONTINUE 

SET  UP  READER  FOR  LIST 
LR=LRDROV (LIST) 

C SET  FLAGS 

LE VEL=0 
ICOUNT= 1 
IEND=0 
IW0RD=0 
IG0=0 
JEND= 1 
10  CONTINUE 
IGO= IGO+ 1 
JG0=O 

C ADVANCE  THROUGH  LIST 
X=  ADVSWR ( LR , K ) 

IF  ( K ) 100,20,100 
20  IF  (LEVEL-LCNTR(LR) ) 150,30,7C 
30  IF  ( N AMTSTt  X ) ) 60,40,60 
40  IF  ( L ISTMT ( X ) ) 50,10,50 
C WE  HIT  A SUBLIST 
50  CONTINUE 

LEVEL=LEVEL+1 
IEND=Q 
IWOR  D=Q 
GO  TO  10 

C WE  HIT  A DATUM  ELEMENT 
6j  CONTINUE 
IENC=0 

IF  (IG0.LE.3)  GO  TO  10 
IF  (JGO.EQ.l)  IC0UNT=IC0UNT+1 
IP  (JGO.EQ.l)  GO  TO  10 
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IWORD= I WORD+ 1 
C STORE  DATUM 

AWORD ( I WORD ) =X 
GO  TO  10 

C WE  HIT  AN  END  OF  SUBLIST 
70  CONTINUE 

LEVEL=LEVEL-1 
IEND= I END+1 

IF  ( IENDoLT.2)  GO  TO  75 
JG0  = 1 
GO  TO  20 
75  CONTINUE 

IF  ( IWORD.GTol)  GO  TO  80 
IF  ( IWORD.LE.O)  GO  TO  20 
C STORE  NUMERICAL  COEFFICIENT  (CONSTANT  TERM) 
SO  LTN ( ICGUNT)=-AW0RD(1) 

IW0RD=0 
GO  TO  20 
80  CONTINUE 

C STORE  MATRIX  COEFFICIENT 
DO  90  1=1,78 

IF  ( AWJRD(2).NE.  ZEE ( I ) ) GO  TC  90 
AM  AT ( ICOUNT, I ) = AWORD ( 1 ) 

IW0RD=0 

GO  TO  (20,160),  JEND 
90  CONTINUE 
95  CONTINUE 

CALL  SLPERR ( 10H  BUCKSV  ) 

100  IF  (LEVEL-LCNTR(LR)  ) 150,120,110 
110  CONTINUE 

LEVEL=LEVEL-1 
GO  To  100 
120  CONTINUE 

CALL  RCELL(LR) 

150  CONTINUE 
JEND=2 
GO  TO  80 
160  CONTINUE 

INVERT  MATRIX  TO  FIND  NUMERICAL  SOLUTION 

CALL  MAT  I NV ( AM AT , ICOUNT, SOL TN , 24, 1 , DET ) 
IF  (DET.EQ.O.)  GO  TO  95 
RETURN 
END 
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SUBROUTINE  COMBY ( N » I » IOUT ) 


C 

C THIS  ROUTINE  CALCULATES  THE  NUMBER  OF  COMBINATIONS  OF 
C N ITEMS  TAKEN  I AT  A TIME 
C WRITTEN  BY  MA  HIRSCHBERG 
C DECEMBER  1975 
C 

IOUT=FACTRL(N)/(FACTRL(N-I )*FACTRL(I) ) 

RETURN 

END 
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SUBROUTINE  EVALP ( AP , FFRM , SOL, NOM,NTOT, I VAL ) 


THIS  ROUTINE  EVALUATES  PI  TERMS 
WRITTEN  BY  MA  HlRSCHBERG 
NOVEMBER  1975 


COMMON/M AXM IN /X AM AX, XAMIN 

DIMENSION  AP(6, 100 ), FFRM ( 100 ) , SOL (24) , I VAL (24) 
DIMENSION  EXPN0(24) , INDX( 24) 

DIMENSION  ARAY( 24) 

C 

C AP  ARRAY  WITH  PRIMITIVE  NAME  AND  UP  TC  5 VALUES 

C FFRM  ARRAY  WITH  NAME  OF  SOLUTIONS 

C SCL  ARRAY  WITH  NUMERIC  SOLUTIONS 

C NFRMS  NUMBER  OF  FORMULAS 

C IVAR  NUMBER  OF  VARIABLES 

C IVAL  NUMBER  OF  VALUES  FOR  EACH  VARIABLE 

C 

C IF  AP  HAS  AN  ENTRY  WITH  THE  VALUE  -9999,9  END  EVALUATION 
C 

WRITE  (6,10) 

10  FORMAT  ( 27H0 VALUE ( S ) OF  PI  TERM  FOLLOW) 

DO  20  1=1, NOM 
EXPNG ( I ) =SQL ( I ) 

20  CONTINUE 
NU  S E = 6 
I 1 =NUSE 
I2=NUSE 
I3=NUSE 
I4=NUSE 
15 =NUSE 
I6=NUSE 
I7=NUSE 
I8=NUSE 
19  =NUSE 
1 10  = NUS  E 
II 1=NUSE 
I 1 2=NUSE 
I 1 3=NUSE 
1 1 4=NUS  E 
I 1 5=NUSE 
I16=NUSE 
II 7=NUSE 
I 1 8=NUSE 
II 9=NUSE 
I ?0=NUSE 
1 2 1 = NUSE 
I 22=NUSE 
I23=NUSE 
1 2 4 = NUS  E 
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ND 1=NUSE 
ND2=NUSE 
ND3=NUSE 
ND4=NUSE 
ND5=NUSE 
ND6=NUSE 
ND7=NUSE 
ND8=NUSE 
ND9=NUSE 
ND10=NUSE 
ND 11=NUS  E 
ND 12  = NUS  E 
ND 13  = NUS  E 
ND 14=NUS  E 
ND 15  = NUS  E 
ND 16  = NUS  E 
ND 17  = NUS  E 
ND 18  = NUS  E 
ND 19  = NUS  E 
ND20  = NUS  E 
ND21  = NUS  E 
ND22  = NUS  E 
ND23  = NUS  E 
ND24=NUS  E 
TERM1=1. 
TERM2=1. 
TERM3=1. 
TERM4=1. 
TERM5=1. 
TERM6=1. 
TERM7=1. 
TERM8=1. 
TERM9= 1 . 
TERM  10  = 1 • 
TERM 11=1. 
TERM  1 2= 1 o 
TERM13=1. 
TERM14=1 • 
TERM15=1* 
TERM16=1. 
TERM17=1. 
TERM18=1 • 
TF  RM 1 9 = 1 • 
TERM20  = 1 • 
TERM2 1= 1 • 
TERM22=1. 
TERM23=lo 
TERM24= 1 • 
IU  SE  = 0 


DO  50  J=1,N0M 
FR  = FFRM ( J ) 

DO  AO  K= 1 , NTGT 
AUSE=AP( 1,K ) 

IF  (FR.NE.AUSE)  GO  TO  40 
IUSE=IUSE+I 
INDX  ( IUS  E ) =K 
GO  TO  50 
40  CONTINUE 

WRITE  (6,42)  FR 

42  FORMAT  ( I5H0NAME  NOT  FOUND, 5X, AID) 

CALL  EXIT 
50  CONTINUE 
KUSE=25-N0M 
DO  150  K = 1 , NOM 
LL  =NGM-K  + 1 
IUSE= INDX ( LL ) 

1 1 = I V AL ( IUSE ) 

IF  ( II.EQ.O)  GO  TO  130 

GO  TO  (124,123,122,121,120,119,118,117,116,115,114,113,112,111, 
1 110,109,108,107,106,105,104,103,102,101),  KUSE 
1C1  11  = 1 1 
NO  1=  1 
GO  TO  130 
1C2  12  = 1 1 
ND2=  1 
GO  TO  130 
1C 3 13  = 1 1 
ND3  = 1 
GO  TO  130 

104  14  = 1 1 
ND4=  1 

GO  TO  130 

105  1 5 = 1 1 
ND5  = 1 

GO  TO  13u 
1C  6 1 6 = 1 1 
ND6=  1 
GO  TO  130 
107  17  = 1 1 
ND7  = 1 
GO  TO  130 
1C8  18  = 1 1 
ND8=  1 
GO  TO  130 
1C  9 19  = 1 1 
ND9  = 1 
GO  TO  130 
llu  110=11 
ND  10  = 1 
GO  TO  130 
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111  111=11 
ND 11=  1 

GO  TO  130 

112  112=11 
NO  12=  1 

GO  TO  130 

113  113=11 
NO  13=  1 

GO  TO  130 

114  114=11 
NO  14= 1 

GO  TO  130 

115  115=11 
NO  15= 1 

GO  TO  130 

116  116=11 
NO  16=  1 

GO  TO  130 

117  117=11 
N017=l 

GO  TO  130 

118  118=11 
ND 13= 1 

GO  TO  130 

119  119=11 
NO  19= 1 

GO  TO  130 

120  120=11 
ND20= 1 

GO  TO  130 
121  121=11 
ND21= 1 
GO  TO  130 
122  122=11 
ND22=1 
GO  TO  130 

123  123=11 
ND23=  1 

GO  TO  130 

124  124=11 
ND24=  1 

130  CONTINUE 

KUSE=KUSE+1 
150  CONTINUE 

0LDTRM=-9999.9 
DO  250  LL  = 1»  1 
DO  248  Jl«ltU»NDl 
IF  (Il.EQ.NUSE)  GO  TO  240 
CALL  NUMBIKK* INDXtli IUSEtKUSE) 
ARAY( 1)=AP(J1+1»IUSE) 
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TERM1=ARAY( 1 ) **EXPNO ( KUSE ) 

DO  248  J 2»  1 » 1 2 » ND2 

IF  ( 12. EQ.NUSE)  GO  TO  240 

CALL  NUMB(KK»INDX,2» IUSE»KUSE) 

ARAY(2)=AP(J2+1,IUSE) 

TERM2=ARAY(2)**EXPN0(KUSE) 

DO  248  J3  = l » 1 3 » ND3 

IP  ( 13. EQ.NUSE)  GO  TO  240 

CALL  NUMB ( KK » INDX » 3, IUSE, KUSE) 

ARAY(3)=AP(J3+1»IUSE) 

TERM3=ARAY(3)**EXPN0(KUSE) 

DO  248  J4  = l » 1 4 » ND4 

IF  ( I4.EQ.NUSE)  GO  TO  240 

CALL  NUMB ( KK » I NCX » 4» IUSE»KUSE ) 

ARAY(4)=AP(J4+l»lUSE) 

TERM4=ARAY(4)**EXPN0(KUSE) 

DO  248  J 5 = 1 » I 5 » ND  5 
IF  ( I5.EQ.NUSE)  GO  TO  240 
CALL  NUMBIKK, INDX, 5, IUSE, KUSE) 
ARAY(5)=AP( J5+1, IUSE) 
TERM5=ARAY( 5 )**EXPNO(KUSE) 

CO  248  J 6s 1 » 16 , ND6 
IF  ( I6.EQ.NUSE)  GO  TO  240 
CALL  NUMBIKK, INDX.6, IUSE, KUSE) 
ARAY( 6 ) = AP ( J6  + 1,  IUSE ) 
TERM6=ARAY( 6 )**EXPNO (KUSE) 

DO  248  J7=l, I7,ND7 
IF  ( I7.EQ.NUSE)  GO  TO  240 
CALL  NUMBIKK, INCX,7, IUSE, KUSE ) 
ARAY(7)=AP( J7+1,IUSE) 
TERM7=ARAY( 7 ) **EXPNO (KUSE ) 

DO  248  J 8 = 1 * I8,ND8 
IF  ( I8.EQ.NUSE)  GO  TO  240 
CALL  NUMBIKK, INDX,8, IUSE, KUSE ) 
ARAY( 8 ) = AP ( J8+1, IUSE) 
TERM8=ARAY(8)**EXPN0(KUSE) 

DO  248  J9=l, I9,ND9 
IF  ( 19. EQ.NUSE)  GO  TO  240 
CALL  NUMBIKK, INCX.9,  IUSE, KUSE) 
AR  AY ( 9)  = AP(J9  + lt IUSE ) 
TERM9=ARAY(9)**EXPN0(KUSE ) 

DO  248  J10  = l,  IlOtNDIO 
IF  ( 1 10» EQ.NUSE ) GO  TO  240 
CALL  NUMBIKK, INDX, 10, IUSE, KUSE) 
AR  AY ( 10)=AP( J10  + 1,IUSE) 

TE  RM 1G  = AR AY ( 10 ) ** EXPNO ( KUSE ) 

DO  248  Jll=l, Ill.NDll 
IF  ( 1 11. EQ.NUSE ) GO  TO  240 
CALL  NUMBIKK, INDX, 11, IUSE, KUSE) 
AR  AY ( 1 1 ) = AP ( Jll  + 1, IUSE) 

TE  RM 1 1 = AR AY ( 1 1 )** EXPNO ( KUSE ) 
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DO  248  J12=l, I12,ND12 
IF  ( 112. EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK, INDX, 12, IUSE, KUSE) 
ARAY( 12)=AP( J12+1, IUSE) 
TERM12=ARAY ( 12 ) **EXPNO ( KUSE ) 

DO  248  J 13*1 » I 13,ND13 
IF  ( 113. EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK,  INDX, 13, IUSE, KUSE) 
AR AY ( 13)=AP( J13+1, IUSE) 
TERM13=ARAY ( 13 ) ** EXPNO ( KUSE ) 

DO  248  J 14*1 , 1 14, ND14 
IF  ( 114. EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK, INDX, 14, IUSE, KUSE) 
AR  AY ( 14)=AP(J14+1,IUSE) 
TERM14=ARAY( 14)** EXPNO (KUSE) 

DO  248  J 15=1, 115, ND15 
IF  ( 115. EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK, INDX, 15, IUSE, KUSE) 
AR  AY ( 15)=AP( J15  + 1, IUSE) 
TERM15=ARAY ( 15 ) ** EXPNO (KUSE ) 

DO  248  J 16=1 , I16,ND16 
IF  ( 116. EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK, INDX, 16, IUSE, KUSE ) 
AR  AY ( 16) =AP( J16+1, IUSE) 

TE  RM 16  = ARAY ( 16 )** EXPNO ( KUSE ) 

DO  248  J 17=1, 117, ND17 
IF  ( 117. EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK, INDX, 17, IUSE, KUSE) 
ARAY( 17) =AP( J17+1, IUSE) 
TFRM17=AR AY ( 17 )** EXPNO ( KUSE ) 

DO  248  J 18  = 1 , 1 18, ND18 
IP  ( 118. EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK, INDX, 18, IUSE, KUSE) 
AR  AY ( 18)=AP( J18+1.IUSE) 
TERM18=ARAY ( 18 )** EXPNO ( KU SE ) 

DO  248  J 19nl , 1 19, ND19 
IF  ( 119. EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK, INDX, 19, IUSE, KUSE) 
AR  AY ( 19 ) = AP ( J19  + 1, IUSE) 
TERM19=ARAY( 19)**EXPN0(KUSE) 

DO  248  J 20  = 1 , 1 20 , ND20 
IF  ( 120. EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK, INDX, 20, IUSE, KUSE) 
ARAY( 20) =AP( J 20+1, IUSE) 
TERM20=ARAY( 20)** EXPNO (KUSE) 

DO  248  J21=l, I21,ND21 

IF  ( 121. EQ.NUSE)  GO  TO  240 

CALL  NUMB(KK, INDX, 21, IUSE, KUSE) 

ARAY(21)=AP(J21+1»IUSE) 

TERM21=ARAY(21)**EXPN0(KUSE) 
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DO  248  J22'=1»I22» ND2 2 
IP  ( I22.EQ.NUSE)  GO  TO  240 
CALL  NUMB(KK,INDX,22,IUSE,KUSE) 

ARAY(22)=AP( J22+1,IUSE) 

TERM22=ARAY( 22 ) **EXPNO (KUSE ) 

DO  248  J 23=1 1 1 23 1 ND23 

IF  ( I23.EQ.NUSE)  GO  TO  240 

CALL  NUMB(KK,INDX,23,IUSE,KUSE) 

ARAY(23)=AP( J23+1,IUSE) 

TERM23=ARAY( 23 )**EXPNO(KUSE ) 

DO  248  J 24* 1 , 1 24, ND24 

IF  ( I24.EQ.NUSE)  GO  TO  240 

CALL  NUMBtKK, INDX,24, IUSE,KUSE) 

ARAY(24)=AP( J24+1,IUSE) 

TE  RM24=ARAY ( 24)**EXPN0(KUSE) 

240  CONTINUE 

TE  RMO  = TERM  1* TER M2*TERM3*TERM4*TERM5*TERM6*TERM7*TERM8 

1 *TERM9*TERM10*TERM11*TERM12*TERM13*TERM14*TERM15*TERM16 

2 *TERM17*TERM18*TERM19*TERM20*TERM21*TERM22*TERM23*TERM24 
IF  (OLDTRM.EQ.TERMO)  GO  TC  246 

WRITE  (6,241)  (ARAY(III),III-1. NOM ) 

241  FORMAT  ( 34H0VALUES  USED  IN  EVALUATING  PI  TERM/ 

1 ( 1H  ,6E20.8 ) ) 

RECIP=1./TERM0 

WRITE  (6,245)  TERMO,RECIP 

245  FORMAT  ( 18H0EVALUATED  TERM  = , E20. 8, 2X, 13HRECI PROCAL  = ,E20.8 

246  CONTINUE 

OL  DTRM=TERMQ 
248  CONTINUE 
250  CONTINUE 
RETURN 
END 
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SUBROUTING  NUMB ( KK» INDX»  LOOP  » IUSEt KUSE ) 


WRITTEN  BY  MA  HIRSCHBERG 
NOVEMBER  1975 

DIMENSION  INCX(!24) 

KUSE=LOOP 

IUS  E= I NDXdNUS  E ) 

RETURN 

END 
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SUBROUTINE  ONEMR ( NSOLS , IRAYS , NCT, I XTRA, I BYPS ) 


THIS  ROUTINE  iPBRMUTES  THE  INPUT  EQUATIONS 
WRITTEN  BY  MA  HIRSCHBERG 
DECEMBER  1975 


DIMENSION  NSOLS ( 720 ) ? IRAYS ( 20 ) 

DIMENSION  I USE  (<20  ) 

DATA  IOK/O/ 

C 

C NSOLS  ARRAY  SHOWING  UNIQUENESS  OF  PERMUTATION 

C IRAYS  ARRAY  WITH  NEXT  PERMUTATION 

C NCT  NUMBER  OF  PERMUTATIONS 

C I XTRA  NUMBER  OF  ITEMS  IN  ARRAY  IRAYS 

C I BYPS  TOTAL  NUMBER  OF  PBRMUTATIONS 

C 


IF  (NCT+WGT.720)  NCT*-1 
IF  (NCT+l.GT.IBYPS)  NCT=-1 
IF  (NCT.ECU-1)  RETURN 
IT  IMES=0 


XTRA=  I XT R*A 


10  CONTINUE 

ITIMES=ITIMES+1 
IF  ( ITIME9.LE.1000)  GO  TO  15 
NCT=-1 
RETURN 
15  CONTINUE 

DO  20  1=  lt-IXTRA 
IlDSE ( I ) = I 


2G  CONTINUE 

DO  50  1=1 »IXTRA 
30  CONTINUE 

XUSE=URAN31( I0K)*XTRA+1. 

IF  ( XUSE.G'T.  XTRA  ) XUSE?XTRA 
IPUT=XUSE 

IF  ( IUSE(TPUTWLE.O)  GO  TO  30 
IRAYS  ( I)  = LPUT 
IUSE  ( I PL|T  1^0 
50  CONTINUE 
IT  EST=0 

DO  60  1 = 1 *. IXTRA 
ITEST-10*ITEST+IRAYS(  I ) 


60  CONTINUE 

DO  70  1=  1 t-NCT 

IF  ( ITEST.lQ.NSOLSl I ) ) GO  TO  1C 
70  CONTINUE 
NCT=NCT+ 1 
NS'OLS  (NCT  I=,IT EST 


C 


RETURN 

END 
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SUBROUTINE  PRMTE( ITI ME , NTOT, SOL ) 


THIS  ROUTINE  MANIPULATES  PI  SOLUTIONS 
WRITTEN  BY  MA  HIRSCHBERG 
NOVEMBER  1975 

COMMON/M AXM I N/X AM AX, XAMIN 
C0MM0N/S0LSAV/S0L1 
DIMENSION  SOL ( 25) 

DIMENSION  SOL  1 ( 25  ) 

GO  TO  (10,30,50,70,90,110),  ITIME 

SQUARE  NUMBERS 
10  CONTINUE 

DO  20  1= 1 , NTCT 
S0L1(  I ) =S0L ( I) 

SOU  1)=2.*S0L(  I ) 

2C  CONTINUE 
RETURN 

CUBE  NUMBERS 
30  CONTINUE 

DO  40  1=1, NTOT 
SOU  I ) =3o  *SQL  1 ( I ) 

Ayj  CONTINUE 
RETURN 

SQUARE  ROOT 
50  CONTINUE 

DO  60  1=1, NTOT 
SOL ( I)=.5*S0L1(  I) 

60  CONTINUE 
RETURN 

CUEE  ROOT 
70  CONTINUE 
A=l./3. 

DO  80  1=1, NTCT 
SOU  I ) =A*S0L  1 ( I ) 

80  CONTINUE 
RETURN 

FIND  MAX  AND  MIN 
90  CONTINUE 

XAMAX=-10000. 

XAMIN=10000. 

DO  95  1=1, NTOT 
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IF  (SOLim.EQ.O.O)  GO  TO  95 
XAMAX* AM AX1 ( XAMAX ,ABSF ( SOLI ( I ) ) ) 
XAMIN=AMIN1(XAMIN,ABSF(S0L1(I)) ) 
95  CONTINUE 

DO  100  1*1 » NTOT 
SOL(  I )=S0L1(  D/XAMIN 
100  CONTINUE 
RETURN 
110  CONTINUE 

DO  120  1=1, NTOT 
SOLI  I )=S0L1( I ) /XAMAX 
120  CONTINUE 
RETURN 
END 
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SUBROUTINE  T ABLUK ( NTOT , SOL TN ) 


WRITTEN  BY  MA  HIRSCHBERG 
NOVEMBER  1975 

THIS  PROGRAM  COMPARES  PI  THEOREM  SOLUTIONS  WITH  LAND  NUMBERS 

DIMENSION  SOLTN ( 24 ) 

DIMENSION  SOL ( 24 ) 

DIMENSION  I TAB ( 93 ) 

DIMENSION  ATABL ( 5 16 ) 

DATA  ITAB/  2, 4, 10, 2, 2,1, 3, li, 10, 13, 2, 1,5, 1,2, 2, 1,2, 1,1, 1,1,1, 

1 4,28,18,6,3,6,4,3,1,1,2,2,1,1,2,1,1,1, 

11, 1,1,1, 1,1, 1,1, 1,1, 1,1, 5, 17, 2, 3, 2, 1,2, 1,3, 1,2, 1,2, 1,1, 1,1, 1,1, 

1 6 , 8, 1,2, 1,1,  1,1,  1,1,  7,4,1, 1,1,1,  8y3ylylyl/ 

DATA  (ATABLII) , 1=1,1331/ 

* -1., 1. ,6HCR0CC0,5HDEBYE,8HEINSTEIN, 

*7HKNUDSEN»5HLAVAL» 4HMACH , 4HNAZE,6HSARRAU,1QH  SMO  LUCKOWS , 

*7H  SP  HEAT,  -2. , 2. , 6HCAUCHY , 5HHOOKE , -1. ,2. ,4HEVAP,4HELAS,  l.yl.y 
*4HHALL,-1. ,-l. , 1. ,9HARRHENIUS,4HEVAP,5HJAC0B,7HPRANDTL,5HRGSBY, 
*7HRUSSELL,7HSCHMIDT»  8HSURF  VEL , 6HV I SOELA S , 7HC0LBURN, -1. ,1. ,1., 
*4HBI0T, 10HBODENSTEIN,8HREYNOLDS,3HCAP,3HNN1,7HNUSSELT,6HPECLET, 

*6H  PCS NOV ,7HPRANDTL»8HSHERW00D,8HSTRQUHAL*  7HTH0MS0N ,5HTRUNC , -o  5 , 

*-. 5, lo , 10HB0USSINESQ,6HFR0UDE,-2. ,2. ,2. ,3HCAP,-2. ,-l.  , 1 . , 5HEKMAN, 
♦4HELAS, 5HEULER,7HFANNING,3HNN2,-2. ,2.  ,4, , 6HTAYL0R , -1.  ,-l.  ,2.  , 
*6HECKERT,6HFR0UDE,-3. ,-l. , 1. , 4HFL0W , 10HMASS  RAT  10,1. ,1. ,1. , 

*9H M PRANDTL.-2. , 1., 1. , 5HREECH , 7HF0UR I ER , 5, .5, 1.  ,7HPRANDTL,-lo  , 
*1.  , 4j , 8HR AD  PRES, -.3333333333, a 3333333333, 1. , 5HSACHS , -2 . , 1 . , 3. , 
*7HGAL ILEO,-. 66 6 6666666,. 3333333333, .3333333333, 10HKIRP I TCHEF, 
♦-1..-1., 1., 1. ,6HBANSEN,7HBINGHAM,7HB0UGUER,7HNUSSELT,8HHEDSTR0M, 
*9HKIRPICHEV,9HK0SS0VICH,8HLAGRANGE, 1CHCR I SPATION, 9HD AMKCHLER , 

*4H EL AS, 5 HELL  IS, 6HGRAETZ, 8HM0MENTUM , 7HNUSSELT , 8H PARTI  CL E ,4HPLAS/ 
DATA  (ATABLl I ), 1=134,310) / 

*8HPIPELINE,-1. ,-l. , 1. ,2. , 9HBR I NKMANN, 9HDAMK0HLER , 10HP0 I SEU I LLE , 
*10HPUMERANTSE, 10HPREDV0D I TL , 1CHSGMMERFELD , -1 . ,-l. ,-l.  , 1 . , 6HHERSEY, 
*8HELSASSER,7HSTANT0N,-1. ,1. ,1. , 1. , 9HDAMK0HLER , 5 HLEWI S , 6HMERKEL  , 
♦8HREYN0LDS, 10HRICHARDS0N,7HSEMEN0V,-1. ,1. ,1. ,2. ,4HBGND, 5HWEBER, 
*6HE0TVQS,9HMAG  INTER , - 1. , . 5 , . 5 , 1. , 6HALFVEN,6HKARMAN,8HMAG  MAC A, 

5, .5, 1. , 1. ,8H HARTMANN, -2. ,-2. ,-l. ,3. ,3HACC,-3. ,-l. , 1. ,4. , 
♦6HM0RT0N , 7HCAP-B0 Y»-3. ,-2. ,-l. , 1. , 8HL AGRANGE , 10HHEAT  TR ANS , -2. , 1 . , 
*1.  ,2. ,8HHEDSTR0M,-1. ,-.5,-. 5, 1. , 7HC0WL I NG , -. 5 , . 5 , . 5 , 1 . , 7HG0UCHER, 
*8HDERYAGIN,-1. 5, -.5,. 5 , 5HEKM AN , -2. , - 1. , 1. , 1 . ,4HEL AS , -. 4 ,-. 2 , . 2 , 
*1. , 5HEXPL0.-2. , 1. , 1. , l.,8HSURATMAN,-1.5,-.5,.5, 1. .7HSLUSH  T,-.75, 
*-. 75, .5, 1. , 8HSP  SPEED, -3. ,-l. ,-l. ,1. ,4HT0MS,-1. ,-.5 ,. 5, 1. , 
*8HLEVERETT,-2.,-l.,l.,2.,8HMAG  PRES , 1. , 1 . , 1 . , 1 . , 10HM  REYNOLDS, 

*~2. »-2. , -1. , 1. , 6HNEWT0N, -c 6666666666, .333333 33 33,. 66666 666 66, 1. , 
♦7HNUSSEL  T ,-.5, -.5, -.5, 1. , 9H0HNE SGRGE , -5.  ,-3.  ,-l. ,1. ,5HP0WER,-1. , 
*.5,1., 1. ,6HREGIER,-1. ,-l. ,-l. ,1. ,1. , 9HKI RP I CHE V , 8HREYN0 LDS , -1. , 
*-l. , 1. , 1. , 1. ,7HH0DGS0N,7HGRAVITY,3HNN5,-2. ,-l. , 1. , 1. , 1. ,5H JOULE/ 
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on  no  on 


DATA  (ATABK  I) , 1 = 311,516)  / 

*5HDARCY,-2.,-l.,l.,l.,2.,7HSQUEEZE,-2.,-l. ,-la, la ,1. ,6HST0KES, 
*7HF0URIER,-1.,-1.,-1.,1. ,2. ,9HMARANG0NI ,-l. ,-l. ,1.,1. ,2., 

♦7HMAG  DYN,7HJEFFREY,9HDAMK0HLER,-.5,1. , 1, , 1. , 1. 5 ,9HLUNDQUI ST,-2. , 
*-lclo  ,lo  ,3,  ,6HKARMAN,10HARCHlMEDES,-l.,-la,-la  ,.5,1.  ,3HCAP, 

*-l .,-.666666 666,-. 333333333,. 666666666, 1* ,8HJ  FACTOR* 6HC0NDEN, 
*-l., l.,l.,2. ,2a,10HCENTRIFUGE,-l.,-l. » 1. , l, , 3. , 8HCLAUSI US » 

*-l.t -.66666 6666*. 333333333*. 666666666, 1. »8HJ  FACTOR, -1.  ,-.5,l.,lo , 
*1.5,4HDEAN,-l.,-.5,.5, 1. , 1. , 5HFRUEH,-1, , 1, , 1. , 1. , 1. ,6HPECLET, 

1 a , — 1 a , — 1* » 1*  » 1« *2* , 7HBAGNCLD  »-3# » — 1 • * ”1 • ,1*  ,1a *1#  *6HTHRING* 
*9HB0LTZMANN,-la,-l.,-l.,l.,l.,4.,6HSTEFAN,-la,-l.,-.5,.5,l.,2., 
*10HELECTR0VIS,-1. , -1 . , 1. , 1, , 2. , 2. ,9HMAG  FORCE, -1* ,-l. ,-.5,l.,l., 
*1#  * 3HNN3  , - 1#  »-  lo  , — a 5 , — *5, 1#  ,1#  ,3HNN^,  — 2#  ,1*  ,1*  ,1a  ,2a  ,3*  , 7 HGR  AS  HO  F , 

10HC0NDENSATI » — 1 * ,1a, la, la, 

♦ 9HDAMK0HL  ER , -2, ,-2. , -1 . , -1. , 1. , 2. , 2. , 6H0CVI RK,-1.  ,-l. ,-l.  ,1.,1., 
*1. ,2.,8HBU0YANCY,-1. ,- 1. , 1. , 1« , 1. , 1„ ,2a ,3# , 8HRAYLEIGH ,-3# ,-2a,-lo, 
*-la, 1., 1 ., l.,4,,7HMCADAMS,-l. . 5 5 5 , . 5 , 1 . ,2. ,2.  ,8HLYK0UDIS/ 

SAVE  SOLUTIONS 
IG0  = 1 

DO  10  I = 1 , NTOT 
SOL ( I ) =S0LTN ( I ) 

10  CONTINUE 


SORT  ARRAY 
15  CONTINUE 
1 = 2 

18  CONTINUE 

IF  (SOLI I)*GE*S0L(I-1) ) GO  TO  20 
X=  SOL ( I- 1 ) 

SOU  1-1)  =SOL(  I ) 

SnL ( I ) =X 
GO  TO  15 
20  CONTINUE 

IF  ( I.EQ.NTOT)  GO  TO  25 
1 = 1 + 1 
GO  TO  18 


FIND  TABLED  VALUES  FOR  SORTED  EXPONENTS 
25  CONTINUE 
C DO  LuOK-UP 


30 


1=1 


<=1 

CONTINUE 
IF  (I.GT.93) 

IU  SE= I T AB ( I ) 

IF  ( IUS5.EQ.NT0T) 


GO  TO  300 


GO  TO  50 


1=  1 + 1 

JMUL= IT  AB ( I ) 
K=  IUSE*JMUL+K 
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DO  40  11=1, JMUL 
1*1  + 1 

K=K+ 1 TAB ( I) 

40  CONTINUE 
1=1  + 1 
GO  TO  30 

WE  HAVE  THE  RIGHT  NUMBER  OF  ELEMENTS 
50  CONTINUE 
KT  EMP=K 
ITEMP=I+1 
IDX= IT  AB ( ITEMP) 

DO  200  11=1, IDX 
ITEMP=ITEMP+1 
LUSE=ITAB I ITEMP ) 

DO  100  J J * 1 , IUSE 
KUSE=KTEMP+JJ-1 

IF  (ABSF(SOL(JJ)-ATABL(KUSE) ) «LE« • 0001 ) GO  TO  100 
GO  TO  150  ' 
iOu  CONTINUE 

WRITE  (6,110) 

110  FORMAT  (16H0LAND  CANDIDATES) 

LI =KUSE+ 1 
L2  =L 1 + LUS  E“ 1 

WRITE  (6,120)  (ATABL(IJ),IJ=L1,L2) 

120  FORMAT  ( lh  ,BA12) 

GO  TO  210 
150  CONTINUE 

KTEMP=KTEMP+IUSE+LUSE 
200  CONTINUE 
210  CONTINUE 

INVERT  SOLUTION  AND  TRY  AGAIN 
GO  TO  (260, 3C0),  IGO 
260  CONTINUE 
I G0  = 2 

DO  270  11=1, NTOT 
SOL ( 1 1 )— S0LTN(  II  ) 

270  CONTINUE 
GO  TO  15 
SCO  CONTINUE 
RETURN 
END 
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APPENDIX  II 

SAMPLE  PROGRAM  INPUTS 
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NUMERICAL  INPUTS 


VARIABLE 

ENGLISH  UNITS 

METRIC  UNITS 

Force  (F) 

104  Lbs 

4 

444820  x 10  Dynes 

Velocity  (V) 

103  Ft/Sec 

30480  — 
sec 

Area  (A) 

2 2 

10  Ft 

92903.04  cm2 

Density  (6) 

.00237  Lb  ^e- 
Ft 

.0012144  D/ne  Sec 
cm4 

Viscosity  (y) 

„ ,rt-7  Lb  sec 

4 x 10  9 

Ft 

1.915201  Dyne  Sec 

2 

cm 

PROGRAM  INPUTS 


8 

FORCED 

VISCOSITY 

AREA 

DENSITY 

VELOCITY 

FORCE 

LENGTH 

TIME 


444820. 

1.915201 

92903.04 

.00122144 

30480. 

-9999.9 

-9999.9 

-9999.9 


04-9999.9 

04-9999.9 

-9999.9 

-9999.9 

-9999.9 


5 

3 

FORCED=FORCE$$ 

VISCOSITY=FORCE*TIME/LENGTH**2$$ 
ARE A=  LENGTH*  * 2 $ $ 

DENSITY-FORCE*TIME**2/LENGTH**4$$ 
VE  LOC ITY=  LENGTH/T IME  $ $ 


NOMOR 


$$ 
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APPENDIX  III 

INPUT  DESCRIPTION  AND  FORMAT 
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CARD 

TYPE 

DATUM 

NUMBER 

NAME  OF 
DATUM 

DATA  ENTRIES 

DATA 

FORMAT 

TYPE 

1 

1 

NPRMS 

Total  number  of  variables  and 
fundamental  units. 

1 

2 

1 

APRIMS 

(1,N) 

Name  of  variable  or  fundamental 
unit  (there  are  NPRM  Type  2 
cards  read) . 

2 

2 

2-6 

APRIMS 
(2-6, N) 

Values  for  fundamental  units 
(used  by  the  evaluation  routine 
EVALP) . -9999.9  signifies  end 

of  data.  Up  to  five  values  are 
allowed  per  variable.  All  com- 
binations of  input  values  are 
calculated. 

3 

1 

NFRMS 

Number  of  equations 

1 

4 

1 

I VARS 

Number  of  fundamental  units 

1 

5 

L FORMS 

Dimensional  equations  (there 
are  NFRMS  equations  which  are 
read  by  the  SYMBOLANG  routine 
INLIST. 

* 

FORMATS 

Type  1 

Format  (15) 

Type  2 

Format  (A10,  2X, 

5E12.8) 

* 

Dimensional  equations  are  read  by  the  SYMBOLANG  routine  INLIST  as  modi- 
fied to  run  on  the  BRLESC  II  Computer) . See  Appendix  II  for  a sample 
of  inputs.  Each  equation  is  terminated  by  a In  addition,  a card 

with  "NOMOR"  (beginning  in  card  column  1)  and  a card  with  "$$"  (also 
beginning  in  column  1)  are  used  to  terminate  the  INLIST  read  and  termi- 
nate input  to  program  BUCKY. 
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APPENDIX  IV 

SAMPLE  PROGRAM  OUTPUTS 
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PI  THEOREM  SOLVER 
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